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We investigate an estimator to measure the primordial trispectrum in equilateral type non- 
Gaussian models such as k-inflation, single field DBI inflation and multi-fleld DBI inflation models 
from Cosmic Microwave Background (CMB) anisotropies. The shape of the trispectrum whose 
amplitude is not constrained by the bispectrum in the context of effective theory of inflation and 
k-inflation is known to admit a separable form of the estimator for CMB anisotropies. We show that 
this shape is 87% correlated with the full quantum trispectrum in single fleld DBI inflation, while 
it is 33% correlated with the one in multi-field DBI inflation when curvature perturbation is origi- 
nated from purely entropic contribution. This suggests that g"^^^ , the amplitude of this particular 
shape, provides a reasonable measure of the non-Gaussianity from the trispectrum in equilateral 
non-Gaussian models. We relate model parameters such as the sound speed, Cs and the transfer 
^«^ ■ coefficient from entropy perturbations to the curvature perturbation, T-jzs with g'^^^ , which enables 

^sj , us to constrain model parameters in these models once g"^^^ is measured in WMAP and Planck. 






INTRODUCTION 



The statistical properties of primordial fluctuations provide crucial information on the physics of the very early 

universe [1-4] (See [5] for a review) . In the simplest single field inflation models where the scalar field has a canonical 

Q kinetic term and quantum fluctuations are generated from the standard Bunch-Davis vacuum, non-Gaussianity of 

I ' the fluctuations is too small to be observed even with future experiments [6-8] . Thus the detection of non- negligible 

Mh departures from Gaussinaity of primordial fluctuations will have a huge impact on the models of early universe. So far, 

most of the studies have focused on the leading order non-Gaussianity measured by the three-point function of Cosmic 

Microwave Background (CMB) anisotropies, i.e. the bispectrum [9-11]. Especially, the optical method of extracting 

the bispectrum from the CMB data has been sufficiently developed [12-18] (for a more general approach, see [19-21]). 

^N I However, future experiments like Planck [22] can also prove the higher order statistics such as the trispectrum [23-25]. 

J^ The trispectrum gives information that cannot be obtained from the bispectrum [26-28]. In addition, it is possible 

\f^ that the trispectrum can be the leading order non-Gaussianity, that is, even if we do not detect the bispectrum, this 

^^ , does not mean that the primordial fluctuations are confirmed to be Gaussian. 

,—1 ' In the local type non-Gaussian models [29-32], the primordial curvature perturbation C, is modeled as 

o ■ c = G + \f^r\ci - iO) + ^9nT\! (1) 

O ! 5 ^5 

^~~^ '■ where Q is a Gaussian variable [33]. In this model, the bispectrum has maximum amplitude for the squeezed 

^ [ configurations in the Fourie space where one of wavenumbers is small compared with others. In these models, the 

• ^ i; trispectrum indeed gives very interesting tests on multi-field infiation models where there are several Gaussian variables 

^\ ' Q using the consistency relation between the amplitudes of the bispectrum and trispectrum [34] . Also the trispectrum 

n3 , can constrain the cubic-order non-linearities in primordial curvature perturbation, g^^"' that cannot be constrained 

■ - - ' by the bispectrum measurements. Estimators to measure the trispectrum in the local-type non-Gaussianity have 

been developed and the kurtosis based estimator [35] have been used to obtain constraints on the amplitude of the 

trispectrum, -7.4 < 5^^710^ < 8.2 at 95% confidence level from the WMAP 5-year data [36]. 

There are another class of non-Gaussian models. A typical example is Dirac-Born-Infled (DBI) inflation [37] whose 
non-Gaussian property was extensively studied by [38-49] (see also [50, 51] for reviews). In this model, like k-inflation 
[52, 53], the inflaton fleld has non-canonical kinetic term and non-linear derivative interactions can give rise to large 
non-Gaussianity of quantum fluctuations. For current observational constraints on DBI inflation see [54-66]. In 
these models, the amplitude of the bispectrum has a peak typically for the equilateral configuration in the Fourie 
space. The shape of the trispectrum is more complicated. For example, for the bispectrum, the equilateral condition 
ki = k2 = ks completely specifies the shape of the bisepctrum, but this is not the case for the trispectrum. The 
shape of the trispectrum has been analyzed in several infiationary models such as single field DBI infiation [67-70] , 
multi-field DBI infiation [71-75] and the models motivated by effective theory of infiation [76, 77]. 
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Regardless of these efforts, since the form of the trispectrum is generally very complicated, estimators for the 
trispectrum in this class of non-Gaussian models have not been implemented yet so far. It was suggested that the 
form of trispectrum given by 

equil 
cv i) kik2kzk^{ ^'+^^\^^+^' f ^^ ' ^ ' 

represents the shape of the trispectrum in equilateral non-Gaussian models very well. Here the trispectrum of the 
curvature perturbation is defined as 

(C(ki)C(k2)C(k3)C(k4)> - (27r)3<53(ki + ka + ks + ]^^,T^{k^MM, ^4), (3) 

where ^^^{k) is the power spectrum given by (((ki)C(k2)) ~ (27r)^(5('^)(ki + k2)fcj^'^7'^(fci). This trispectrum (2) 
appears in DBI inflation as a contribution from the fourth order interacting Hamiltonian (the "contact interaction" ) 
[69, 70]. In the effective theory of inflation, it was shown that the trispectrum of this shape can have the amplitude 
that is not constrained by the bispectrum measurements [76]. In Ref. [70] , it was suggested that this trispectrum 
can be used to construct an estimator because by introducing the integral 1/M" — (l/r(n)) /(,°° i("~^)e~*'^*, this 
function is factorisable (sec Appendix A). Therefore, in this paper, we compare the shapes of trispectra in single field 
and multi-field DBI infiation with Eq. (2) based on a shape correlator introduced by Regan et.al [78] and investigate 
whether the estimator constructed from the trispectrum (2) represents the shapes of trispctrum in these models or 
not. 

This paper is organized as follows. In section II, we review the shape correlator introduced by Regan et al. [78]. 
In section III, we study the overlap between the shape given by Eq. (2) and trispectra in single field and multi-field 
DBI inflation models. In section IV, we give theoretical predictions for g^L** i'^ some concrete theoretical models. 
We conclude in section V. In Appendix A, we present the optimal estimator using Eq. (2) explicitly. In Appendix 
B, we summarise the shape function of the reduced trispectra appeared in general single field k-inflation models and 
give the shape correlations among the representative shapes. In Appendix C, we check the validity of our method to 
relate the amplitude of the estimator to the theoretical predictions using the bispectrum. 

II. THE SHAPE CORRELATOR 

In this section, we review the shape correlator introduced by Regan et al. [78]. 

A. Shape functions 

First we exploit the symmetry of the trispectrum to define the reduced trispectrum as follows [23] . We rewrite the 
definition of the trispectrum as 

(C(ki)C(k2)C(k3)C(k4))c = (27r)-'' j d^K [6{ki + ka - K)5{ks +U+K) {Tdk^M, ^sM; K) 

+rc(k2, ki, kg, k4; K) + rc(ki, k2, k4, kg; K) + rc(k2, ki, k4, kg; K)) 

+ (k2 ^ kg) + (k2 ^ k4)] . (4) 

Then we need to consider only the reduced trispectrum 7^ from one particular arrangement of the vectors, such as 
7$(ki, k2, kg, k4; ki2) with ki2 = ki -I- k2 and form the other contributions by considering permutations. Here, for 
the later convenience, we use the symmetrised reduced trispectrum 

7^''^'"(ki,k2,k3,k4;ki2) = - [7c(ki, k2,k3, k4; ki2) + 7c(k2,ki, k3,k4; ki2) 

+7c(ki,k2,k4,k3;ki2) + 7c(k2, ki,k4, kg; ki2)] , (5) 

and from now on we omit the superscript sym for simplicity. 

The reduced trispectrum is a function of six variables. We can choose them to be (fci, fc2, ^3, ^4, fci2, O4) where 64 
represents the deviation of the quadrilateral from planarity which is specified by the triangle (/ci, fe, ^12)- We find 



that in terms of these variables ku — |ki + k4| is expressed as 

1 

2k 



fcl4 — K]^ + fc4 o;„2 V^l + ^12 ^2)1^4 +^^12 ^3) 



±^^j4kfkf, - {kl + kl^ - klY^^klkl^ C0S2 04 - [kl + kl^ - kir , (6) 

which imphcs that the vahd range of cos 64 is constrained by 

\cos0,\>^-^^f^. (7) 

2fci2K4 

Motivated by the relation between the CMB trispectrum and the trispectrum for (, the shape function for the 
reduced trispectrum is defined as 

Sr{ki,k2,k3,k4,ki2,04) = {kik2k3k4)'^ki2T(:{ki, k2, ks, ki; ki2,04) . (8) 

Regan et al. [78] proposed to define an overlap between two different shape functions Sj- and 5*7-' as 

F{Sr,Sr') = dVk d{cose4)Sriki,k2,k2,k4,ki2,94)Sr'{ki,k2,k3,k4,ki2,94)w{ki,k2,k3,k4,ki2) , (9) 

where w is an appropriate weight function. The weight function should be chosen such that S'^w in k space produces 
the same scaling as the estimator in I space and we adopt the one used in Rcf. [78], 

w{ki,k2,k3,k4,ki2) = - — , — , , w, , , — —j—:- (10) 

kl2[kl + K2 + Kl2j(fc3 + K4 + ki2) 

With this choice of weight, the shape correlator is defined as 

CiSr, Sr^) = F{Sr,Sr>) ^ 

^F{Sr,Sr)F{Sr',Sr') 

B. Parameterisation of six parameters 

First to parameterise the magnitude of the momenta, we use the semiperimeter of the triangle formed by the vectors 
ki,k2,ki +k2, 

q = -{kl + k2 + ki2) ■ (12) 

From the scaling behaviour, the form of the shape function on a constant-g cross section becomes independent of 
q and we can write 

Sr{ki,k2,k3,k4,ki2,64) = /(g)5*r(fci, ^2, A:3, A:4, fci2, 6*4) , (13) 

where ki = ki/q and fci2 = ki2/q. Since we are restricted to the region where the momenta {ki, k2, ^12) and (fcs, ^4, fci2) 
form triangles by momentum conservation, we will reparameterise the allowed region to separate out the overall scale 
q from the behaviour on a constant q cross-sectional slice. This five-dimensional slice is spanned by the remaining 
coordinates. For triangle (fci,fc2,fci2) we have 

fci2 - 9(1-/?), (14) 

kl = |(l + a + /3), (15) 

k2 = |(l-a + /3), (16) 



while for triangle (^3,^4,^12) 



ki2 = €qil-S), (17) 

fc3 = |(l + 7 + '5), (18) 

fc4 = f(l-7 + <5), (19) 



where e parameterises the ratio of the perimeters of the two triangles, i.e. e = {k^ + ^4 + fci2)/(fci + ^2 + fci2)- We do 
not lose the generality to consider 1 < e < 00. The different expressions for ki2 imply that 



l-/3 = e(l-(5), (20) 



from which 5 is eliminated to give 



^^3 - |(-l + /?+(2 + 7)e), (21) 

ki = |(-l + /?+(2-7)e). (22) 

The conditions for triangle {ki, ^2, ^12) that < fci, ^2, fci2 < <Z imply that < (3 < 1 and —(1 — /3) < a < 1 — /3, 
while the condition for triangle (^3,^4,^12) that < ^3,^4,^12 < efci2 imply that —(1 — /3)/e < 7 < (1 — /?)/e- 
Furthermore, in terms of these variables, the condition (7) is expressed as 

|cos.4|>'^^^^-T;^-^ + ^^\-^^IV^^'. (23) 

' '- (l-/?)(-l + /3 + (2-7)e) ^ ' 

In summary, we have the following domains, 

0<g<oo, l<e<oo, 0<^<1, -{l- 13) <a<l- P, — < 1 < ^-^ , (24) 

e e 

together with Eq. (23). In practice, we introduce a cutoff for the integration of e as S-rS-rw is decreasing with e 
asymptotically after integrating out the dependence of a, 7 and /3 in the overlap integral, we set the cut-off to be 
e = 10 but the dependence on this cut-off is very weak. Also we should emphasize that the CMB measurements will 
never prove the parameter region where e ^ 1. 

Making use of this parameterisation, the shape function (13), the weight function (10) and the volume element can 
be rewritten as 

Sr{ki,k2,k3,k4,ki2,04) = /(q)5r(Q:,/3,7, e,6'4), 

w{ki,k2,k3,k4,ki2) = —— -T- , (25) 

4e(l-/3) 

dVk — dkidk2dk^dk4dK — eq dqdad/Sd'jde . (26) 

It is worth noting that although the integration dVk is five-dimensional in Eq. (9), for scale-invariant shape functions 
with constant f{q), it is enough to evaluate shape correlations only for the four dimensional slices with constant q. 

III. SHAPE CORRELATIONS 

In this section, we study the overlap between Eq. (2) and the trispectra in single field and multi-field DBI inflation. 
It is worth mentioning that from the definition of the shape correlator (11), the shape correlations are independent 
of the normalisations of shape functions. We will discuss the normalisation of the trispectrum in the next section. 

A. Equilateral shape 

First, we find that the shape function for the trispectrum (2) is given by 

ge^ml ^ ^eqmlgc,^ (27) 

o 

where S!ji is given by Eq. (B9). We will assume the scale independence of the spectrum P^ in the rest of the paper. 
In k-inflation model, this class of models are characterised by P^4x ^ X~'^P^xx , X~^ P^xxx for the action given 
by Eq. (Bl). It was also shown that, in the context of the effective theory of inflation, it is possible to construct 
consistent inflationary models where the trispectrum is characterised by this shape function and its amplitude is not 
constrained by the bispectrum [76]. 
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FIG. 1; We show the e dependence of S^"* 5^"* w (sohd line), Sj- ("> g'^^'' yj (dashed hne) and Sj- '■''> g^^'' ^ (dot 
ted line) by integrating out the dependence of a, 



and 7 in the overlap integration. We find that S^ S^"' w and 



-iD B I((t) Qcquil 



Sj- y") g'^^'' y; asymptote to oc 1/e'' while S^"' S^"' w oc 1/e^ for large e. We find that while Sj- ("> g'^'^^ yj jg always 



-yDBI(s) Qequil 



positive, Sj- Sjf w become negative above some critical value of e. Because of this, the full overlap between S- 



-lequil 



DBI(s) 
T 



and 



Sj^ takes much smaller value than the one estimated at configurations with e = 1. We have normalised so that the values 
become 1 at e = 1. 



As this shape function depends on a, /?, 7, e, first, we clarify the e dependence of the signal which is given by 
w. We find that after integrating out the dependence of a, j3 and 7 in the overlap integration, the 



-tequil Qequil 



"T 



s: 



T 



amplitude of the signal is proportional to 1/e^ asymptotically. This shows that the dominant contribution to the 



equil Qequil 



W. 



signal for this shape is coming from e ~ 1. In Fig. 1, we show the e dependence of Sj^ Sjf 

Next, we examine the a, /3, 7 dependence of S'^™ . For this purpose, we plot 5^"* (q;,/3,7, e) evaluated at e = 1 
for given /3 in Fig. 2 where S*^"* is symmetric under the exchange of a and 7 for the configurations with e = 1 and 
the physical region is given by 1 — /3>a,7>/3 — 1. 



B. Single filed DBI inflation 

It was suggested in [69] that the shape function corresponding to the reduced trispectrum of single field DBI inflation 
at leading order in the slow-roll expansion is given by 



QDBI{a) 
bj- 

pfDBI{a) 






-3S^ 



1 
64' 



S: 



T 



1 
64 



^ 64 



r 



,6^4 



06 c; 



(29) 

(30) 



where S^\ S^ and S^' are given by Eqs. (B12), (B13) and (B18), respectively. 



Similar to the case of the equilateral shape, first we examine the e dependence. In Fig. 1, we plot 5'^ 



DBI{(t) Qequil 



'r 



w. 



We find that after integrating out the dependence of a, /3 and 7 in the overlap integration, the amplitude of the 
signal is proportional to 1/e'^ asymptotically. This shows that the dominant contribution to the overlap between the 
single field DBI model and the equilateral shape is coming from e ~ 1. The difference of the asymptotic e dependence 
between the trispectra corresponding to the differences of the shapes between S^ , which is coming from the contact 
interaction and S^ , S^ , S^ , which arise from the scalar exchanges. This was pointed out in Ref. [70] by considering 
the double squeezed limit (fca — k^ — ki2 -^ 0). Therefore, it is natural that the asymptotic e dependence between 
S'^"* and Sj- is different, as Sj- is obtained by a linear combination of S'. 

Next, we examine the a, (3, 7 dependence of Sj- 

e — 1 for given (3 in Fig. 3 where Sj- is symmetric under the exchange of a and 7 for the configurations with 

e ~ 1 and the physical region is given byl — /3>a,7>/3 — 1. Except for the region with very small value of /3 



-T- , Oj- and ijj- . 
For this purpose, we plot Sj- (a, /3,7, e) evaluated at 



-^equil 



IS 



(from to ~ 0.1) the shapes are very similar to 5^"* . This explains that the overlap between Sj- and S'. 

sufficiently large for the configurations with e — 1. 

Table I provides a summary of correlations between S'^"* and Sj- . In addition to the correlation considering 

full configurations dealing with five dimensional parameter space, for comparisons, we also consider the configurations 




FIG. 2: In this group of figures, we plot S^'"' (a,/?,7, e) with e = 1. /? = 0,0.1,0.25,0.5,0.75,0.9 for upper left, upper 
right, middle left, middle right, lower left, lower right, respectively. Vertical axes for each plot denote a and 7. For these 
configurations, S^*" is symmetric under the exchange of a and 7 and the physical region is given by 1 ~ /? > q:,7 > /3 — 1. 
We set iV'='J"'' = 1. 



limited with e = 1 and equilateral configurations (e = 1, a = 7 = 0). The overlap decrease if we include the non- 
equilateral configurations keeping e = 1. This is due to the difference of the shapes for small /3. Also the different 
asymptotic behaviours with respect to e further reduces the overlap if we integrate over e. However, even after we 
perform the all integration, the overlap remains high at around 87% level. 




FIG. 3: In this group of figures, we plot S°'^"'''\a,'y, P) with e = 1. /3 = 0,0.1,0.25,0.5,0.75,0.9 for upper left, upper 
right, middle left, middle right, lower left, lower right, respectively. Vertical axes for each plot denote a and 7. For these 
configurations, Sj- is symmetric under the exchange of a and 7 and the physical region is given by 1 — /? > a, 7 > /? — 1. 

We set Ar-Dsi(<^) ^ i 

C. Multi-field DBI inflation 

It was suggested in [73] that the reduced trispectrum of multi-field DBI inflation dominated by originally purely 
entropic perturbations at leading order in the slow-roll expansion is given by 



^DBI(s) 



JSfDBIis) 
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rS': 
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DBI{s) _ 
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,6^4 



(j)'^C 



tT, 



ns ' 



(31) 
(32) 



where S^ and S^ are given by Eqs. (B31) and (B32), respectively. 

Again, we first examine the e dependence of the overlap (in Fig. 1). We find that after integrating out the dependence 



of a, /3 and 7, the amplitude of the overlap, 5'^"* Sj- w is proportional to 1/e^ asymptotically as in the single 
field case. This shows that the dominant contribution to the signal for this shape is coming from e ^ 1. 

The asymptotic e dependence of Sj- ^''' is the same as that of Sj- '^' . We find that asymptotically S*^*", S- 



S2d 

T 



and Sfp'' given by Eqs. (B15), (B17) and (B20) give the dominant contribution to both Sj- and Sj- , which 

characterises the asymptotic e dependence. However, as is shown in Fig. 1, Sj- '^'^ > g^"^" ^ become negative above 
some critical value of e. This reduces the final overlap once we integrate over e. 

Next, we examine the a, /3, 7 dependence of Sj- . For this purpose, again we plot Sj- (a, /3, 7, e) evaluated 

at e = 1 for given /3 in Fig. 4 where Sj- is symmetric under the exchange of a and 7 for the configurations with 

e — I and the physical region is given byl — /3>q;,7>/3— 1. 




FIG. 4: In this group of figures, we plot S^'^"-'\a,^,l3) with e = 1. /3 = 0,0.1,0.25,0.5,0.75,0.9 for upper left, upper 
right, middle left, middle right, lower left, lower right, respectively. Vertical axes for each plot denote a and 7. For these 
configurations, Sj- is symmetric under the exchange of a and 7 and the physical region is given by 1 — /3 > 0,7 > /3 — 1. 

We set N°Bi{s) ^ ^^ 



Except for the region with very small value of /3 (from to ^ 0.1) the shape is different from S'^"* . This explains 
that the overlap between Sj- and S^^^ is not so large even for the configurations with e = 1. 

Table I provides a summary of correlations between S!^^^ and Sj- . As in the single field case, in addition to 
the correlation considering full configurations dealing with five dimensional parameter space, for comparisons, we also 
consider the configurations limited with e — 1 and equilateral configurations {e — 1, a = 7 = 0). 

Table I shows that the overlap becomes smaller once we include the non-equilateral configurations with e — 1. 
This is clear from the shape difference for /3 > 0.1. In addition the shape correlation for full configurations becomes 
further smaller once we integrate over e. As explained before, this is due to the fact that while Sj- ^'^> g^^'^ ^ 
is always positive, Sj- v^J^egui ^ become negative above some critical value of e. This confirms the fact that the 
shape dependence of trispectrum can in principle distinguish multi-field DBI inflation models form single field DBI 
inflation models shown by Refs. [72, 73]. In practice, the overlap still remains at 33% level after integrating over all 
the shape parameters and the equilateral shape could still be used to get a reasonable estimation for the constraints 
on multi-field DBI inflation model. 

IV. THEORETICAL PREDICTIONS FOR gX" 

In this section, making use of the shape correlations investigated in the previous section, we give theoretical 
predictions for g^^^ ■ As a consistency check, we applied the same method to estimate the amplitude of the bispectrum 
in DBI inflation in Appendix C. 

In k-infiation models, as is shown in Appendix B by setting P4X ^ X^^P^xxx,X^^P^xxx, the shape function is 
given by 

rrl2p 

qequil ^^ -* ,4X QCi {'i^\ 



IGP^i^s 



equil 



Then, by comparing Eqs. (27) with (33), we find g^^ is obtained as 

where we have used P^ = H^ /{'\:XcsP^x) for single field k-infiation. 

In order to express the amplitude of trispectrum in single field DBI inflation in terms of g^L** ' ^^ rewrite Eq. (29) 
in the following form: 



equil 
gDBIi.) ^ Q^^ 3 5nL 



r 64 r g4 r 54 r 



(35) 



where the numerical factor in Eq. (35) is chosen so that when g'^^^ = 1 and C{Sj- ^'^ , 5^™ ) = 1; the following 
conditions arc satisfied, 

7-,/^r>B/(CT) QDBI{a)\ jp/ Qequil riequil\ 

F{S^^'^"\S^^'''') > 0. (36) 

Of course, C(Sj- , S*.^"* ) = 1 is not true in reality and this factor will enhance the amplitude of the signal for a 

given 5n'l'* ■ This term is necessary because when we use the estimator related with 5^"* for the signal whose shape 
is characterised by Sj- , the observed signal is suppressed by C{Sj- , S*.^"* ) and it is necessary to compensate 

this. Then, by comparing Eqs. (29) with (36), we can relate g^L** with the sound speed c^ as 



equil 



17 



9nl = -4 ' (37) 



where we have used 7^$ = H'^/{2(jP) for single field DBI inflation. 

Similarly, in order to express the amplitude of trispectrum in multi-fleld DBI inflation in terms of .g^'iT* ' ^^ rewrite 
Eq. (31) in the following form: 



equil 

^DBI(s) _ , 2-p3 ffNL 

' "^ /t"/ c'-C*-B/(s) qequil\ 



CC2 t QSl , QS2 , QS3 

8^^ + 576^^+64^^ + 192^-^ 



(38) 



10 

where the numerical factor in Eq. (38) is again chosen so that when g^L^' — 1 and C{Sj- , S'^™ ) = 1, the fohowing 

conditions are satisfied, 



F{S. 



DBI{s) nDBI{s) 



T 



\s. 



T 



equil Qequil\ 



\ rpfriequti ciequtL\ 



DBI{s) Qequil\ 



^(5;^"^^'^5: 



r 



> 0. 



(39) 



Then, by comparing Eqs. (31) with (39), we can relate 5^^** with the sound speed c^ and the transfer coefficient 
that relate the amplitude of original entropy perturbations to the final curvature perturbation T^g as 



equil 

9nl 



2.2 



(40) 



where we have used P(j = H T^g / {2(j) ) for multi field DBI inflation. 



equil 



In Table I, we summarise theoretical predictions for gj^^ for th^ models discussed in this section 





Overlap- full 


e= 1 


equilateral 


theoretical prediction for g^'j^'' 


requil 

/nl 


equilaterai shape 


1 


1 


1 


{ZX\lP,ix)/{l&P,x) 


requil 

/nl 


single DBI 


0.87 


0.90 


0.92 


17/4 


-0.36/4 


multi DBI 


0.33 


0.60 


0.85 


2.2/(c*Tts) 


-0.36/(c?Tts) 



TABLE I: Shape correlations against 5^"* for full configurations, the configurations restricted to e = 1, equilateral configura- 
tions (e = 1, Q = 7 = 0) in the model with equilateral shape motivated by effective theory of inflation, single field DBI inflation 
and multi-field DBI inflation. We also summarise theoretical predictions for g^'i!'* ^'I'l /nl"* ^^ these models. 

It is instructive to compare the values (37) and (40) with previous results of tnl [73] based on the matching of the 
amplitude at a speciflc equilateral conflguration. We deflne the non-linear parameter 



Tnl 



rc(ki,k2,k3,k4)fe?fc|fcifc|P^-3 



[{kiki + kikD{k. 



13 



-14 



('1,31,3 I ^3jL3Vl. 



12 + ^13 ) 



('^l ^3 + '^2^4)1^12 +^14] 



(41) 



and evaluated it for the conflguration specifled by ki = k2 = k^ = k^ = k, ki2 = fcis = fci4 = {2/\/3)k. For single 
field and multi-field DBI inflation models, we obtain 



Tnl 



TNL 



0.56 
0.12 

^s-^TZS 



for single fleld DBI inflation 
, for multifield DBI infiation 



(42) 
(43) 



Using the same procedure, we get tnl = '^9'nh^ /9a/3- By comparing this to (42) and (43) we can estimate g'^™ as 



■7 4 4 

equil ~c.-L 

5nl - "X 



for single field DBI infiation 



equil 

5nl ~ 



0.94 
^47-2 



for multifield DBI infiation 



(44) 
(45) 



which underestimates the amplitude by factor 2^4 compared with the results of Eqs. (37) and (40). This demon- 
strates that unlifce the bispectrum case where the matching of the amplitude at the equilateral configuration gives a 
reasonable estimation for /^l"* , it is necessary to calculate the overlap between the shapes in full five-dimensional 
parameter space to extract the amplitude of the trispectrum g'^™ ■ 



CONCLUSION 



It is well fcnown that there are many interesting early universe models that predict equilateral type primordial 
non-Gaussianity motivated by string theory and effective field theory. Taking into account the fact that future 
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experiments such as Planck can prove even next order statistics, it is important to study the primordial trispectrum 
in these models. For example, we had shown previously that the trispectrum can in principle distinguish multi-field 
DBI inflation models from single field DBI inflation models fi:om the shape dependence [69, 72, 73]. On the other 
hand, at the practical level, since the form of the trispectrum is too complicated, the estimator in this class of models 
had not been implemented explicitly. 

Therefore, in this work wc have presented a method to estimate primordial trispectrum in equilateral type non- 
Gaussian models such as k- inflation model whose action is given by Eq. (Bl), single field DBI infiation and multi- field 
DBI inflation. Our method is based on the following two facts. One is that the equilateral shape given by Eq. (2) 
becomes factorisable by introducing the integral 1/M" = (l/r(n)) /p°° i(n-i)e~^* as was suggested in Ref. [70]. The 
other is that in terms of the shape correlation proposed by Ref. [78] , we can relate the amplitudes of trispectra with 
different shapes. 

After reviewing the shape correlator, we have calculated the overlaps between the equilateral shape and the shapes 
of trispectra in single fleld DBI inflation and multi-fleld DBI inflation. We have shown that the shape is 87% correlated 
with the one in single field DBI inflation, while it is 33% correlated with that in multi-field DBI inflation when the 
curvature perturbation is originated from purely entropic perturbations during inflation. We have summarised the 
overlaps including the configurations restricted to e = 1 and equilateral configurations (e = 1, a = 7 = 0) in Table I. 
We found that the main difference between the the equilateral shape and the shape in single field DBI inflation comes 
from the conflgurations with /3 ^ 0, which can be seen even in the equilateral conflgurations (fci = k2 = ks = k^). For 
the shape in multi-field DBI inflation, as the behaviour of the shape function is different from the one in single field 
DBI inflation [69, 72, 73], the overlap becomes smaller. Regardless of this, when we take into account of the fact that 
this overlap is calculated in the five-dimensional parameter space, the 33% correlation is not necessarily small. For 
example, the overlap between equilateral shape and local shapes, which depend on cutoffs in the integration due to 
divergences in various limits, is less than 2%. 

Then, we have given theoretical predictions for g"^^^ , which enables us to constrain this type of non-Gaussian 
models from future experiments. For the model with equilateral shape motivated by k-inflation, we obtained g^L^' = 
(3X3c2p4x)/(16Px), while for single field DBI infiation and muhi-ficld DBI inflation, 5^!™' == l^/c* and .g^^™' == 
2.2/(c^r.^^), respectively. To obtain this value, instead of matching the amplitudes of the shape functions at a specific 
point in the parameter space, we have adopted an overlap function, F{S-t-, S-j-'), defined in Eq. (9), which involves 
integration over five-parameters. 

Before closing, let us comment on the detectability of the trispectrum in future experiments. According to the 
estimation in Refs. [16, 76], the observational errors on /nl and ^nl scales as 

Vy N : VcN ■ 

C pix ' s pix 

where iVpix represents the number of data points of the experiment. Therefore, current limit on (/nl by WMAP is of 
order 10 , while the future experiments like Planck [22] and 21-cm line experiments [79] are expected to produce a 
limit (7nl ^ 10^ and (/nl ~ 10^, respectively. 

For the models like single field DBI infiation where non-Gaussian parameters are given by /nl ^ c^"^ and (/nl ^ c-'^'^i 
the trispectrum is not detectable even by the Planck satellite since there is already a constraint like c^ > 10~^ from 
/nl ^ 10^ from the bispectrum measurement [76]. However, for multi- field DBI infiation models where non-Gaussian 
parameters are given by /nl ~ c~'^T^g and (/nl ^ cj'^T^g, it might be possible to detect the trispectrum if there 
is a large transfer from the entropy mode. For example, it is detectable by Planck if T-jis ~ 10. In this context, to 
construct a concrete theoretical model which gives a large transfer coefficient T-jis is important. The model with the 
equilateral shape motivated by effective theory of inflation [76] can give (/nl much larger than 10^ which is detectable 
by Planck while keeping the value of /nl to be just of order one. 

Finally, although we have not studied in this paper, it is known that the ghost inflation also gives equilateral type 
non-Gaussianity [80, 81] and recently the shape dependence of the trispectrum was also calculated [82, 83]. It might 
be interesting to express the amplitude of the trispectrum in terms of the estimator proposed in this paper. 
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Appendix A: Optimal estimator 



In this section, we present the optimal estimator to detect the trispectrum in equilateral type non-Gaussian models. 
As was shown in the main text, Eq. (2) is a representative form of the trispectrum which has sufRciently large overlaps 
between trispectra in physically motivated models such as single field and multi-field DBI infiation. Moreover, this 
trispectrum can be written in a factorisable form as 



T^(fci,fc2,fc3,fc4) 



equil 

_3nl 



fclfc2fc3fc4(^i±H^^^^)5 ^ 

equil ^5 
^""^^ Vl I dttU{t,ki)d{t,k2)d{t,k3)d{t,k4), 



where 



d(k,t) = — exp(—kt). 
k 



(Al) 
(A2) 

(A3) 



Then the connected part of the trispectrum of the CMB temperature anisotropics is calculated as 



{aiirmai^rn2airm3aumi)c = TTj'Pc / d^Yi^rniYi^rniYhmsY^rm t dt r drji^{r,t)"/i^{r,t)"/i^{r,t)ji^{r,t) 



where 



24 



7h{r,t) = - / dkikfgTii{ki)F{ki,t)ji^{ki,r), 



(A4) 



(A5) 



gxii is the radiative transfer function, ji. is the spherical Bessel function and w;. is an experimental window function. 
The optimal estimator is given by [78] 



equil *^ 

5nl ~ ^' 



where 



S ~ -;j7 2_^('^hmiai2m20-l3m3ai^„ii)c (C a)l^rni{C a)l^m2{C a)i^m^{C a)i^ 



lirui 



'0(C )limil2m2\^ O'Jhmsi^-' OJ;4m4 + ^(G )limil2m2[^ Jhmslim^i 



(A6) 

(A7) 
(A8) 



and 



^ ny,\'^liniiO,l2rn2^l3rn30,l4rn4)c\0-l'^'m[^l'^rn'2^l'^rn'^(''l'^m,'^)c\'^ )lirni,l[m[['^ ji2m2;'2™2^ )hm3,l'^m'^['^ )lim4,l'^ 



Here 



Cl-m-^l-mj — {o-limi(^ljm.j) , (C ajj^m. — Ci j ttl- 



Using the expression for the trispectrum (A2), the estimator for the equilateral trispectrum can be written as 

45 equtij,3 r r r 

5 = li|^ [t'^dt IrHr [dn 



where 



h 

(D(ri,r, i)^)MC = y^^Wl,Wl.^l^^l^{C~'^)l,„i^J^rn,Yi.miYi.mj. 
hli 



(A9) 
(AlO) 

(All) 

(A12) 
(A13) 
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The ensemble average can be evaluated using the Monte Carlo simulation. The fisher error bound of g*^^^ is given 
by F~^ where 

where we assume the covariant matrix is diagonal and the reduced trispectrum is given by 

^(3/4 (^) = ''^24 '"^'" / ^^^^ / '^^(^'^'^hWi2WisWu-fi,{r,t)-fi^{r,t)-fi^Xr,t)ji^{r,t). (A15) 



Here h^u, is given by 



Ku. ^ ^jSh±2m±lM:±ll '. i . ,Ai6) 



Appendix B: Shape functions in general single field k-inflation 

Here, based on our previous work [69] , we summarise the shape functions for the reduced trispectra in general single 
field k-inflation described by the following action: 

S=^Jd^x^^[R + 2P{X,<P)], (Bl) 

where (j) is the inflaton field, R is the Ricci scalar and X = — (l/2)(7^'^9^09i/0, where g^jy is the metric tensor. 

For this class of models, the third and the fourth order interaction Hamiltonian of the field perturbation S(j) in the 
flat gauge at leading order in the slow-roll expansion are given by 

Hf^\v) = I d?x AaScf)'^ + BaScj)' {dS<j)f , (B2) 

H^iv) ^ fd^x[l3iS(P'^ + l32S(f>'^dS(Pf+l33{dS(pf^, (B3) 

where prime denotes derivative with respect to conformal time 77 and coefficients A, B, /3i, P2 and P3 are given by 

/o V / 2 \ \/2X 

A = -^^[Pxx + ,:XP^xxx], B = ^^Pxx- (B4) 



2 V 3 

A = Pxx (1 - Icl^ 2XP.XXX (1 - \ci) + ^P'xxx \x'Pax. 
P2 = -\P.XX [l - \ci^ + ^XclPxxx, 

Then, the shape function S!^ is composed of two parts 

qk qk(cont) qk(scalar) IV>.f\\ 

where S'j- " denotes the contribution from the contact interaction and 5^ ""^ denotes that from the scalar 
exchange interaction, respectively. 



iryk(cont) . ■ n 

bj> IS given by 



^^Hcont) ^ (_24/3ic35f^ - /32C,5f^ - 2p,c;\S!^-) ^N» . (B7) 
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Here S^ , S^ and Sj^ arc the following shape functions: 



S^ 



S^ 



3 perms. , 



S!^ 



(Eti k. 

fci2fc^fc|(k3 ■k4) 

fci2fcgfc|(ki ■ ka) 



fci2(ki-k2)(k3-k4) L ^^»<r^ 




1 + 3(^3 + ^4)^^^ fc3/C4 



1 + 3 



(B8) 



(B9) 



3 perms. , (BIO) 



where "3 perms." denotes the permutations (fci o ^2), (^3 ^ fc4) and (/ci O fej^s ^ ^4)- In Eq. (B7), N = 

r,. .1 1 rik(scalar) . . -, 

bimilarly, d^ is given by 



„2 tt2 atIO 
oX 



(Bll) 



Here Sj^ , S^ and S*^ are the following shape functions: 



S'^ 



-9ki2ikik2k3k4y/^ 



Pi{ki,k2, -ki2, fcs, k4, ki2) - Ti{-ki, -fe, -ki2, k^^k^, ^12) 



-J"l(fc3, ^4, -fci2, fci, ^2, fci2) - J'i(-fc3, -fc4, -fci2, ^3, k^, fcia) 



3 perms. . 



(B12) 



5f^ 



s-f^' 



S'f^<= 



S'f^'' 



C'-52a _|_ QS2b _|_ QS2c _|_ QS2d 

= -fci2(fcifc2fc3fc4)'/'(kl • k2)(k3 • k4) 



(B13) 



J"2(-fcl2, fcl, ^2, fcl2, fcs, ^4) - -^2(-fcl2, -fcl, -^2, ^12,^3, k^) 



-J"2(-fci2, fcs, ^4, fci2, fci, fc2) - p2{-ki2, -^3, -^4, fci2, fci, fc2) 



3 perms. 



(B14) 



-2A:i2(fciA:2A:3fc4)'/'(ki • k2)(ki2 • k4 



J"2(-fci2, fci, fc2, ^3, fc4, fci2) - ^2(-fcl2, -fcl, -^2, ^3, /C4, ^12) 



+-^2(^3, ^4, -fcl2, fcl2, fcl, ^2) - -^2(-fc3, -^4, -fcl2, ^12,^1, fc2 



3 perms. 



(B15) 



2ki2{klk2k3kiy/^iki2 ■ k2)(k3 • k4 



J"2(fci, fc2, -fci2, fci2, fcs, fc4) - p2{-kl, -fe, -fcl2, ^12, ^3, ^4) 



-J'2(-A:i2, ^3,^4, fci, ^2, ^12) - -^2(-fcl2, -^3, -^4, ^1,^2, ^12) 



3 perms. 



(B16) 



4fci2(A:ifc2A:3A:4)'/'(ki2 • k2)(ki2 • k4) 



J"2(fci, fc2, -fci2, ^3, /C4, fci2) - ^2(-fcl, "^2, "^12, ^3 , ^4, ^12) 



-^"2(^3, ^4, -fcl2, ki,k2,ki2) - J"2(-fc3, -^4, "^12, ^1,^2, ^12) 



3 perms. 



(B17) 
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S^^ 



jj- -T uj- -^ uj- -^ Uj- , 



5f^» = 3fci2(fcifc2fc3fc4)'/'(k3-k4) 



-^3(^1, ^2, -^12, fcl2, fca, ^4) - -^3(-fcl, -^2, -^12, fcl2, ^3, ^4) 



+J"4(-fci2, fc3, fc4, fcl, ^2, fcl2) - -^4(-fcl2, -^3, -^4, ^1,^2, fc 



12 J 



3 perms. , 



5:^3. ^ 6fci2(A:iA:2fc3fc4)'/'(ki2-k4) 



-^3(^1, ^2, -^12, ^3, ^4, fci2) - Tai-ki, -k2, -ki2, ks, ki, fci2) 



+^"4(^3, fc4, -ki2, fci, ^2, fci2) - Pi(-k^, -ki, -ki2, ki,k2, fci 



3 perms. , 



^:^3e ^ 3fci2(A:ifefc3fc4)'/'(ki-k2 



J^i{-ki2,ki,k2, k3,k4, ki2) - J^A{-ki2, -fci, -fc2, ^3, kA, ki2) 



+J"3(/c3, ^4, -A:i2, fci2, fci, ^2) - fsi-ks, -ki, -ki2, ki2, fci, ^2) 



3 perms. , 



S^^^ = -6fci2(fcifc2fc3A:4)'/'(ki2-k2) 



J^i{ki,k2, k^i2,k3, ki, ki2) - J"4(-fci, -fe, fc-12, ^3, ^4, fcl2) 



+J'3(A:3, ^4, -^12, ki,k2, ki2) - Tsi-ka, -k^, -ki2, ki,k2, k 



12 J 



3 perms. 



(B18) 



(B19) 



(B20) 



(B21) 



(B22) 



where again "3 perms." denotes the permutations (ki ^ ^2), (^3 "H- k^) and (ki 4->- A:2,/c3 ^ /c4). Here we have 
defined four Ti functions (with i = 1, . . . , 4) as follows; 



1 1 / A A^ 

J^l{ki,k2,k3,k4,k5,kG) = -4|fcifc2fc3fc4fc5fc6h .3^3 ( 1 + 3-T + 6— y 



(B23) 



J^2{ki , k2, k^, k4, k^, ke 



\k1ki\2 1 
k2k^k^k(i\^ AC 



fcs + fee ^k^ke 



A 



A" 



+ ^ (k2 + ks + ^5 + fee + -^ ((^2 + ^3) (fc5 + fee) + 2fc5fc6) + 2:^^^^fctM 
2 / 1 

+ ^(A:5fc6 + (fc2 + fc3) (^5 + ke) + k2k3 + -j (fc2fc3 (fcg + /jg) + 2/25^6 (^2 + A:3)) 

k2k-ik^k^^ 6 



-24 



k2kzk^kQ 



) + ^ (feA:3 (fc5 + fee) + ^5^6 (fc2 + ^3) + 2^-^^^^^ 



T3{ki,k2,k3,k4,,k5,ke) = 2^ 
F4,{ki,k2,kz,ki,k^,kQ) = 2^ 



C4 

fcifc2fc3fc4|5 1 

Ifcsfcel" -4C3 
kik/^k^kQ\^ 1 

|fc2fc3|^ ^'C 



k5 + kfi ^k^ka 3 /, , ^k^ka 
1 + . + 2^n^ + - I fcs + fee + 2 



yt 



^^2 C 



yt 



12- 



.^5^6 



(B24) 



(B25) 



A A^ fc2 + fc3 ^ ^(fc2 + fc3) + fc2fc3 



_4 

+ 3^ (^ (fe + fc3) + 2fc2fc3) + 12/C2fc3 



(B26) 



where ^ is defined by the sum of the last three arguments of the Ti functions as ^ = fc4 + k^ + fcg and C is defined 
by the sum of all the arguments as C = fci + ^2 + ^3 + ^4 + ^5 + ^6 • 

In Tables II, we summarise the correlations among these shape functions for the full configurations, the configura- 
tions with e = 1, the equilateral configurations (e = 1, a = 7 = 0), respectively. 



It is worth noting that the following properties 

F{S!^, Sr) = J2a,F{S!ji, S'r) ^ 

i 

F{St, St) = 5]a,a,F(^V, 5f ) . 



(B27) 
(B28) 



«j 



16 





Cl 


C2 


C3 


Sl 


S2 


S3 


§2 


S3 


Overlap- full 


1.00 


-0.55 


0.35 


0.98 


0.83 


-0.89 


-0.056 


0.24 


e = 1 


1.00 


-0.64 


0.72 


1.00 


0.86 


-0.92 


0.32 


0.17 


equilateral 


1.00 


-0.72 


0.79 


0.98 


0.90 


-0.93 


0.48 


0.44 



TABLE II: Shape correlations between the equilateral shape and the shapes of primordial trispectra in general single field 
k inflation. The correlations between the shapes that appear in multi-field inflation models are also shown in the last two 
columns. 



hold for the shape function given by 



Sj- — 2_^ aiSj 



T y 



(B29) 



where i ~ ci, C2, C3, si, S2, S3 and a^'s are corresponding coefRcients. By combining the shape correlations obtained in 
Table II and the properties shown above, we can calculate correlations of the shape functions in any general single 
field k-inflation models against the equilateral shape (B9), once the action (Bl) is specified. 

EspeciaUy, in the case of single field DBI inflation, the coefficients in the Hamiltonians (B2) and (B3) are given by 



A^-- 



1 



/3i 



1 

2c7(/)2 



B = 



1 



/32 



20cf 
1 
4c? (/>2 



/33 



8c, 



In multi- field DBI infiation model [73] , in addition to the shape functions S'Si , 
convenient to define the following shape functions S^ and S^ given by 



(B30) 



S^, S^, S^' 54?, 54?, we find it 



^T- 



'r- 



'T ■ 



'r ■ 



s-^ 



5f^" 



S-^" 



S-^' 



S-^" 



QS3 QS3a QSSb 



QSSc QSSd 



(B31) 
(B32) 



The table II shows the shape correlations between these shapes and the equilateral shape. It is clear that these 
shapes have significantly low correlations which explain the reason why the final correlation between the shapes in 
the equilateral model and the multi-field DBI models is lower than the single field model. 

Appendix C: Bispectrum estimation for single field DBI inflation 

In this section, we explain our method to compare the amplitudes of the bispectrum using the overlap integration. 
We introduce the shape function Sb 



SB{ki,k2,h) = hhh^ciki, fc2, ^3) , 
for the bispectrum of the curvature perturbation i3^(fci, fc2, k^) defined by 

(C(ki)C(k2)C(k3)) = {27TfS(''\k^+k2+ks)Bc{k,,k2,k3) 



(Cl) 



(C2) 



In [43], it was shown that at leading order in the slow roll expansion and small sound speed limit, shape function 
for the bispectrum in the single field DBI infiation is given by 



S 



DBI 



\ i>3 



2i^2 E ''^'''3 



-Ykf 



(C3) 



where K ~ ki + k2 + k^. 

The shape (C3) is not factorisablc and it is not easy to perform an optimal analysis using CMB observations. Thus 
a factorisablc shape which approximates (C3) was proposed by [13] which is given by 



nequil ^2 ^ requiljj t 2 ( 



1 



1 



1 



1 



pL3 



7,3 p 



1,3)1,3 



1,2 L2 1,2 
"'l"'2''-3 



"'I "^9 rCr) 



5 perms 



^2 "-3 



(C4) 
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where the permutations act only on the last term in parentheses. 

One way to relate /^l"' to the prediction of DI 
the equilateral configuration ^1=^2 = ^3. We get 



One way to relate /^l"' f o the prediction of DBI inflation is to compare the amplitude of the shape function for 



/Nr'=-Y^- -0.32c-. (C5) 

In fact, in [20], these two shapes are shown to be 99% correlated to each other based on the following primordial 
shape correlator for two different shape functions Sb and Sb' 

r fQ Q \ Fb{Sb,Sb') .„ . 

\/Fb{Sb, Sb)Fb{Sb',Sb') 

which is constructed from 

Fb{Sb,Sb') ^ dUkSB{ki,k2,k3)SB'{ki,k2,h)wB{ki,k2,h) . (C7) 



In Eq. (C7) the integration is performed for the region where the triangle condition for (fci, fc2, k^) holds and weight 
function wb is given by 

WB = , ^ ,\ , . (C8) 

ki + k2 + K3 

For these two slightly different shapes, it would be enough to match the amplitudes of the bispectra evaluated at the 
equilateral configuration (ki = k2 = fca) where the amplitude have a peak by setting the relation (C5). However, it is 
more appropriate to match the amplitude and shape of the bispectra by taking into account all possible configurations. 
In this context, we use a different way to estimate the amplitude of the bispectrum in single field DBI inflation using 
the information of all possible configurations. For this purpose, we rewrite Eq. (C3) in the following form: 



requil 1 / 1 1 1 \ 



CB{S'^''',ST'')'^^k^ \ K^^ ' ^ 2if2 



where the numerical factor in Eq. (C9) is chosen so that when /^l"* — 1 ^nd Cb{S§^^ , S'^^^ ) — 1, the following 
conditions are satisfied, 

771 (ciDBI cDBI\ 771 / riequil nequil\ 

I'Bi^B i'-'B J - ^b(^b ^^B )y 

FB(5f^^5B''"') > 0. (CIO) 

Of course, Cb{S^^^ , S'^"^ ) = 1 is not true in reality and this factor will enhance the amplitude of the signal for a 
given /nl"* ■ Then, by comparing Eqs. (C3) with (C9), we can relate /nI"* with c^ as 

/X' = -0.36c;2 , (Cll) 

where we have used V(^ = H'^/{2(p^) for single field DBI inflation. As is expected, this gives almost the same value as 
the one given by Eq. (Cll), due to the fact that there is a large overlap between the two shapes. However, for the 
trispectrum the difference between the two approaches tends to be large as the trispectrum has five parameters even 
assuming the scale invariance and thus matching the amplitude at a specific point in the five-dimensional parameter 
space is generally not enough to ensure that we get the same signal. In this case, it is more appropriate to use all the 
shape information using the overlap integration to compare the shape and amplitude of trispectra. 
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